#R commands Intro SNA - Filip Agneessens

getwd()
setwd(“E:/workshop”)
setwd("E:/workshop")
getwd()
list.files()
a<-5.264
a
A
b=3.2
b
a*3
a+10
c<-a+10
c
ls()
c==a
c!=13
d<-(c!=13)
d
rm(a)
ls()
rm(a,b,c,d)
ls()
?rm

install.packages("sna")
library(sna)
mat1<-matrix(c(0,0,0,0,1,0,1,1,0,0,0,1,0,0,1,0),4,4,byrow=TRUE)
mat1
NROW(mat1)
n<-NROW(mat1) 
gplot(mat1)
gplot(mat1, displaylabels=TRUE, arrowhead.cex=2, vertex.cex=3)
sum(mat1)
gden(mat1)
sum(mat1)/(n*(n-1))
rowSums(mat1)
colSums(mat1)
degree(mat1, gmode="digraph", diag=FALSE, cmode="outdegree", rescale=FALSE)

mat1[2,1]
mat1[2,1:3]
mat1[2,]
t(mat1)
mat1t<-t(mat1)
mat1+mat1t
sum(mat1+mat1t==2)
sum(mat1+mat1t==1)
sum(mat1+mat1t==0)
dyad.census(mat1)
grecip(mat1, measure = c("edgewise"))
grecip(mat1, measure = c("dyadic.nonnull"))
grecip(mat1, measure = c("dyadic"))
mat1*mat1t
mat1%*%mat1
mat2<-mat1
mat2[1,4]<-1
mat2[2,4]<-0
mat2

gplot(mat2, displaylabels=TRUE, arrowhead.cex=2)
mat2%*%mat2
mat2%*%mat2%*%mat2
triad.census(mat2)
geodist(mat2)
GM<-geodist(mat2)
GM$gdist
GM$gdist<Inf
M.reach<-(GM$gdist<Inf)
M.reach
M.reach*1
M.gdist<-GM$gdist
M.gdist[M.gdist %in% c(Inf)]<-NA
M.gdist
M.gdist[M.gdist %in% c(NA)]<-9
M.gdist

setwd("E:/workshop")
list.files()
data1<-read.table("klas03e-friends-waveA.dat")
data1
gden(data1)
degree(data1, gmode="digraph", diag=FALSE, cmode="outdegree", rescale=FALSE)
outd1<-degree(data1, gmode="digraph", diag=FALSE, cmode="outdegree", rescale=FALSE)
hist(outd1)
table(outd1)
ind1<-degree(data1, gmode="digraph", diag=FALSE, cmode="indegree", rescale=FALSE)
hist(ind1)
table(ind1)

dyad.census(data1)
grecip(data1, measure = c("edgewise"))
triad.census(data1)
geodist(data1)
list.files()
attr<-read.table("klas03e-demographics.dat")
attr
attr[,1]
gender<-attr[,1]
gplot(data1, displaylabels=TRUE,arrowhead.cex=2,vertex.cex=3,vertex.col=gender)